Asymptotic stability and periodic solutions of 
vector Lienard equations 
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Abstract 



We prove the asymptotic stability of the equilibrium solution of a class 
. of vector Lienard equations by means of LaSalle invariance principle. The 

key hypothesis consists in assuming that the intersections of the manifolds 
in {V = 0} be isolated. We deduce an existence theorem for periodic 
solutions of periodically perturbed vector Lienard equations. 
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1 Introduction 

In this paper we are concerned with the vector Lienard equation 

g : X" + f(X)-X' + g(X) = 0, (I) 

where / e C°(Q, ET), open connected subset of 1", g 6 Lip([tf =1 (a h h), Et n ), 
the space of locally lipschitzian functions on n™=i( a «i^); an d ■ denotes the in- 
ner product of IR™. We assume the origin O to be an isolated equilibrium point 
, of the equivalent system 

X' = Y, Y' = -g(X) - f(X) ■ Y, (2) 

In the planar case (n = 1), the equation (fT]) has been widely studied from many 
points of view: stability of equilibria, boundedness and oscillation of solutions, 
existence, uniqueness or multiplicity of periodic solutions (see [3] , [S] , [5] ) . The 
higher dimensional case present substantial additional difficulties, which prevent 
straightforward extensions of planar results. This is in particular true for an 
equilibrium's stability properties, which are essential in studying the dynamics 
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of perturbed sytems. Even if such properties are often analyzed by linearizing 
the system at the equilibrium, this cannot be done if some component of f(X) 
or some partial derivative of g(X) vanish at such a critical point. 

Among the few available results, we point out two papers. The first one is 
concerned with the study of the singular perturbation problem ex + R(x)x + 
g(x) = for < £ < 1, performed in 0]. For small positive values of £ the 
solutions approach the subspace x = —R (x)g(x). 

A class of Lienard vector equations was considered in pQ, under quite re- 
strictive hypotheses. In [T] g(X) is assumed to be the gradient of a suitable 
positive definite homogeneous scalar function, f(X) to be the jacobian matrix 
of a homogeneous vector field F(X) on H™, such that F(X) ■ g(X) be positive 
definite. This ensures that the conditions required by Barbashin-Krasokskii- 
LaSalle theorem ([6], thm. 79) hold for the system 

X' = Y-F(X), Y' = -g(X), (3) 

which is equivalent to (HJ. In fact, such hypotheses are satisifed by taking the 
energy V(X, Y) — \Y\ 2 /2 + G(X) as a Liapunov function, since the derivative V 
of V along the solutions of ^ vanishes only on the set {X = 0}, i. e. on the "Y 
axis" . The non-positive invariance of such a set is an immediate consequence of 
X' = Y on the on the set {X = 0}. 

In this paper we study the asymptotic stability of the origin for the system 
[[2"]) under different hypotheses, admitting the possibility of a larger vanishing 
set {V = 0}, possibily containing several hypersurfaces not contained in the set 
{X = 0}. This allows to apply our result to equations excluded by the approach 
of [I], as 

f x 'i =2/1 f Vi = -xi~yi{xi~ x 2 ) 2 ,^ 

1 x' 2 = y 2 \ y' 2 = -x 2 -y 2 (x 1 +xl) 2 . 

We assume g(X) = (g\ (xi), . . . ,g n (x n )), so that the coupling of the n scalar 
equations in (fT]) is entirely due to the dissipative terms f(X) ■ X'. Then the 
system ([2]) appears as a Lienard-like perturbation of a natural Hamiltonian 
system with separable variables. Actually, the above Liapunov function V is the 
hamiltonian function of the unperturbed system. We prove the equilibrium's 
asymptotic stability assuming that the intersection of the manifolds in V — 
consists of isolated points (theorem [TJ . In order to do this, we prove that every 
solutions starting at a point of any of such manifolds leaves it immediately, 
proving the non-invariance of V = 0. In corollary[T]we extend our result to prove 
the equilibrium's global asymptotic stability. In corollary [21 we apply theorem 
fT] in order to prove the existence of periodic solutions for small perturbations of 
the type 

X" + f(X)-X'+g{X) = h(t,X,X',s), (5) 
where h(t, X, X', e) = h(t + T, X, X', e) e C* 1 (1R 2ti+1 x [0,e),]R™), T > 0. 
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2 Lienard equation, n > 1 



Let us set X = (xi,...,x n ), Y = (yi, . . . , y n ), Z = (xi, . . . , x n , yi ...,,y n ). Let 
us consider a map / = (/i, ...,/„) denned on an open connected subset 57 of 
IR™ containing the n-origin Ox- For a scalar map g : (a, 6) — » IR, we say that 
,g G Lip((a, b), IR) if g is locally lipschitzian in all of (a, &). 
Let us consider a system of n coupled Lienard equations 



x'l + fi(xx,...,x n )x' 1 +gi(xi) =0 

%2 + M^l, • • • , XjijaJj + 32(^2) = 
X, n ~\- f n (*^1 j • * * ) ^n)^-n "i - .9^ (^n) 0: 



(6) 



where / € C°(Q,]R), ^(a;,) G Lip((a,i, h), H), Oj < < &i, i = 1, . . . , n. The 
equivalent first-order system is 



yi 



y n 



( „,i 



y[ =-9i(xi)-yih{xi,...,Xn) 



(7) 



We assume <?i(0) = 0, Xigi(xi) > for a;, G (a i; 0) U (0,&i), i = l,...,n, so 
that the 2n-origin O = (Ox,Oy) is an isolated equilibrium point of 0. Also 
in this case, if the c/i's are differentiable at 0, in order to study the stability of 
O one could compute the eigenvalues of the linearized system at O, which is 
non-degenerate if and only if n"=i5i(0) 0. Moreover, the eigenvalues' real 
parts depend on the values of fi(Ox), i — similarly to what happens 

for n = 1, 

,1,2 _ -fi(Ox)±y/fi(Ox) 2 -4gm 

Xi - 2 ' 

Hence, if fiiOx) vanishes and 3^(0) > for some 1 < i < n, a different approach 
is required. This occurs also if <7-(0) = 0, for some 1 < i < n. 

The Liapunov function approach is more general, in two ways. It applies to 
degenerate critical points, and allows also to give an estimate of the equilibrium's 
region of attraction. In order to apply LaSalle invariance principle to a suitable 
Liapunov function V, we have to prove that the vanishing set {V = 0} is not 
positively invariant. Such a set is the set-theoretical union of hypersurfaces. In 
next theorem we prove that, if such hypersurfaces' intersections are isolated, 
then {V — 0} is not positively invariant. From now on we use the word isolated 
referring to the topology of IR". For definitions abopu dynamical systems, we 
refer to [2]. 

Theorem 1 If 

1) Xigi(xi) > for x t G (a,, h), Xi ^ 0, i = 1, . . . ,n; 

2) fi(X)>0, in 0, i = l,...,n; 
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3) for every subset S C {1, ...,n} the set {xi = 0, for i e S} fl {/j 
I), for i ^ S} consists of isolated points; 

4) ri™=i{/i — 0} consists of isolated points; 
then O is asymptotically stable. 

Proof. Let us set 

V(xi,...,x n ,yi,...,y n ) = ^2 [ G i( x i) + y) > ( 8 ) 

where Gi(x) — J^ 1 gi(s)ds, with i = 1, . . . ,n. By the above assumptions, V is 
positive definite at O. In order to study the stability of the equilibrium point 
O, we study the sign of the derivative of V along the solutions of (0) 

n 
i=l 

The hypothesis 2) implies that V(Z) < in a neighbourhood of O, hence the 
origin is stable. In order to apply the LaSalle principle, we study the orbits' be- 
haviour on the set W = {J2ti Vih{ X ) = 0}- The equality J^ti Vih{ X ) = 
holds if and only if one of the following holds: 



a) V\ = • • • = Vn = 0, 

b) for some subset S C {1, . 

c) fi = ... = /„ = 0. 



. , n} one has J/, — for i £ S, and fi = for 



Let us denote by Wq, Wq, Wq, the set defined by the respective condition, so 
that Wo = WgUW&UWg. We show that every orbit 7(2) such that 7(0) € W 
leaves immediately Wq, Wq, Wq, hence it leaves immediately Wo- In other 
words, we prove that if 7(0) £ Wo, then for all t £ (0, e) one has 7(4) £ Wo, 
for some positive e. We also set Wq = {yi = for i £ 5, /j = for i £ S}. 

Case a). Let 7(4) be an orbit of 0, such that 7(0) G W§. Since t/i(0) = 0, 
for i = one has y((0) = -gi(xi(0)) in 0. If Zi(0) ^ 0, then t/((0) = 
— gi(xi(0)) 7^ 0, hence 7 leaves immediately the set Wq. As a consequence, 
xi(0) = 0, so that 2/1(0) = 0, and at 7(0) the system appears as follows: 



x\ 



= 
= 

= 



y'i 
y' n 







-52(0:2) 



(9) 



We may apply to j/j, i > 1, the argument just used on j/i, proving that either 
7 leaves immediately the set Wq, or 2^(0) = 0. As a consequence, one has 
2/i(0) = ... = 2/,,(0) = &i(0) = ... = x n (0) = 0. Hence the unique orbit that 
does not leave the set yi = y% = is the equilibrium point O. 
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Case b). Assume 7(0) G Wq. As in case a), we may prove that, if yi(0) = 
for i G S, then Xi(0) — 0, for i G S, otherwise 7 leaves the set {yi = 0, for i G S}, 
hence it leaves immediately the set Wq. If 7 remains in Wq, then it remains in 
the set {xi = 0, for i G S} D {/i = 0, for i ^ S 1 }. By the hypothesis 3), such a 
set consists of isolated points, hence there exist X G H™, such that A(t) = X, 
for i > 0. As a consequence, one has yi(t) — x'^i) = 0, i = l,...,n, hence 
Wq C Wq. Now we may use the fact that ^(t) leaves immediately Wq, hence 
also Wq. 

Case c). Let 7(0) G W Q . By hypothesis 4), n" =1 {/i = 0} consists of isolated 
points. If j(t) starts at a point of Wq and remains in Wq, then its x-components 
are constant. Then one has yi{t) = x'^t) = 0, for i = 1, . . . , n, hence on Wq the 
system has the following form 

A = f y[ = -gi(xi) 

(10) 

<4 = I Vn = -9n{x n ). 

Now we can repeat the argument of case a), proving that y[(0) = . . . = y' n {0) = 
0, hence, at 7(0), Xi(0) = ?/i(0) = 0, i = 1, . . . , n. This shows that 7(0) remains 
in W Q only if 7(0) = 0. * 

We give a couple of examples of systems satisfying the hypotheses of theorem 
[TJ The system 



x[ =y x f y[ = -xi -y x x\{x 2 - \) 2 {x x + l) 2 

x'2 =y-2 \ y'2 = ~ x 2 - y2x\(x x - i) 2 (x 2 + 1) 2 . 

One has {/1 = 0} n {x 2 = 0} = {(0, 0), (-1, 0)} and {/ 2 = 0} n {x t = 0} 
{(0,0), (0,-1)}. 

Moreover {h - 0} n {/ 2 - 0} - {(0,0), (0, -1), (1, 1), (-1,0), (-1, -1)}. 



(11) 



The system 



Vx f V'x = -x 1 -y 1 (xl + 2x 2 .-l) 2 

yi \y'i = -x 2 -y 2 {2x 2 1 +xl-l) 2 . 



(12) 



One has {/1 = 0} n {x 2 = 0} = {(0, 1), (0. -1)} and {f 2 = 0} n {^i = 0} = 
{(1,0), (-1, 0)}. Moreover {h = 0} n {/ 2 = 0} = {(±^, ±^)}. 

The theorem [1] may be extended to prove the global asymptotic stability of 
an equilibrium point. 

Corollary 1 If Q — TR n , (a,, &,) = R for i = 1, . . . , n, and the hypotheses of 
theorem^ hold, then O is globally asymptotically stable for Ffjj. 

Proof. Under the above hypotheses, the Liapunov function V of the theorem[T]is 
dcfmicd on all of HI", with V < on all of IR™. By LaSalle invariance principle, 
every solution 7 tends to the largest invariant set E contained in V = 0. As 
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proved in theorem [TJ E = {O}. This gives the global attractivity of the origin, 
hence its global asymptotic stability. Jft 
As a consequence of theorem [TJ we prove an existence result for periodic 
solutions of small periodic perturbations of the equation [TJ 

Corollary 2 If the hypotheses of theorem\T\ hold, and h(t,X,X',s) S C 1 (H x 
x IR™ x [0,e),IR n ), e > 0, then there exists e* G (0,e) such that the equation 
f5j) has a periodic solution <f>(t,s) for all e £ (0, £*), such that 4>(t,e) tends to 
the null solution as e — > 0. 

Proof. The asymptotic stability of a critical point of an autonomous system is 
equivalent to its uniform asymptotic stability, hence we may apply theorem 15.9 
in [7] in order to prove the statement. Jit 
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